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Abstract 

We develop a Wigner function representation of the quantum transport theory of the conduction 
band electrons in Rashba effect resonant tunneling structures with a phonon bath. In narrow band 
gap heterostructures, spin splitting occurs mainly as a result of inversion asymmetry in the spatial 
dependence of the potential or as a result of external electric field. This "zero magnetic field spin 
splitting" is due to the Rashba term in the effective mass Hamiltonian. The quantum transport 
equations are derived using multi-band non-equilibrium Green's function formulation in generalized 
Kadanoff-Baym ansatz. 
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Zero magnetic field spin splitting |lfl of the conduction band states in bulk and compound 
semiconductors has attracted great interest recently. The possibility of a spin polarized 
current source based on conventional nonmagnetic semiconductors has led many researchers 
to investigate the possibility of obtaining efficient spin-polarized current source, a necessity 
in spintronics. Bulk inversion asymmetry (as in zincblende semiconductors) dominates in 
large band gap semiconductors whereas the structural inversion asymmetry ( the Rashba 
effect ) becomes important in narrow band gap semiconductors. 

The Rashba effect predicts a spin splitting of the conduction band linear in k\\. For 
example, the spin splitting in the lowest conduction subband exceeds 0.02 eV for in-plane 



wave vector k\\ = 0.05 A 1 







A similar splitting occurs in the light hole band proportional 
to fcy and in the heavy hole band proportional to fc| 3|. The case of conduction band 
electrons is more attractive for the spin transport devices since these have a much longer 
spin relaxation time. 

Wigner function modeling of charge transport in multi-barrier resonant tunneling struc- 
tures has been quite popular in the literature due to its success in dealing with the dissipation 
and the open boundary conditions 4]. Similarly, it is expected that one should be able to 
model the spin transport in the tunneling structures using the Wigner function. 

We derive two conduction spin subband Wigner function equations for zero magnetic field 
resonant tunneling structures. These equations are derived for the first time in the literature. 
When the spin-orbit interaction is included in the single conduction band effective mass 
Hamiltonian, the Wigner function becomes a 2 x 2 matrix in spin subspace. The diagonal 
elements fa and fn are the spin up, spin down electron distribution functions respectively. 
The off-diagonal elements fr± and fa describe the spin coherence. Note that fa = fa. 

The Rashba effect resonant tunneling structures have been discussed from both theoretical 



and experimental points of view in the literature. It was shown that an asymmetric dou 
barrier resonant tunneling diode (DB RTD) might provide spin polarization above 50 % |5j 
Koga et al proposed a triple barrier resonant tunneling diode (TB RTD) that can reveal a 
high degree of polarization 99.9 % g|. Recently Ting and Cartoixa proposed the resonant 
inter-band tunneling spin filter |2|. 

If the z- axis is chosen as the growth direction, the effective mass Hamiltonian with the 
spin-orbit term can be written as 



1 Pn a 



where a denotes the Pauli spin matrices, p is the momentum operator, py = (p x ,p y ,0) is 
the in-plane momentum vector. We assume that the in-plane momentum, pn, is conserved 
across the device, a is the spin-orbit coupling constant (Rashba constant) and it is material 
dependent (inversely proportional to the energy gap and the effective mass) and proportional 
to the interface electric field along the growth direction. Therefore it is possible to tune this 
parameter by applying an external bias. The experimental value of the Rashba constant is 
in the order of 10~ 12 eV.m. 

The Hamiltonian for Rashba effect resonant tunneling structure can be written as, 



H = 





H n 




p 2 

Pz 2 J {z) Pz+ 2 J {Z) +E c (z) i a h p_ 

p 2 


(1) 








[ l tP+ P^J(z)P^+ 2m*(*) + E c(z) 





where p z = —ifrj^ , P± = p x ± ip y , and E c (z) = E c + V(z). E c is the conduction band edge 
and V(z) is the self-consistent potential. 
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Throughout the paper, a denotes the spin index and there is a summation over repeated 
indices. The four-dimensional, (3+1), crystal momentum and its conjugate variable lattice 
coordinate are represented as p = (p,E), r = (R, T). Under the assumption that the 
self-energies are slowly varying with respect to the center of mass coordinates, generalized 
Kadanoff-Baym equation for in phase-space-energy-time domain can be written as 



d 1 
ih—G<,(p,r) = —- - I dp 2 dr 2 K c H (p,r - r 2 ;r,p - p 2 )G<,(p 2 ,r 2 



+E>„(p,r)G<,(p,r) - ^(prfGfyfar) (2) 



where 



K%(p,r- r 2 ; r,p-p 2 ) = K\{p,r - r 2 ; r,p-p 2 ) - K A (p,r - r 2 ; r,p-p 2 ) 



and 



dpidriejq)(-pi.(r - r 2 )) exp(-n.(p -p 2 ))A(p ± y,r =f — 



The self-energy function for the electron-phonon scattering in phase-space-energy-time 
representation can be written as 

E>J(p,r) = ±J dqG>J(p + q )D>><( q ). (3) 

Assuming the phonon bath is in equilibrium, the Fourier transforms of the phonon Green's 
functions can be written as, 
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D<(q,E') = 



-ihM$[(N q + 1)5(E' - OJ + N^iE' + fi q )], 



(4) 



D>(q, E) = -ihMlKN^ + 1)5(E' + fi q ) + N q 5(E' - fi q )] 



(5) 



where M q is the electron-phonon scattering matrix element. Therefore, inclusion of the 
phonon scattering gives the following scattering functions 



crcr 



E jjs I dciG^ip + ^E + rjQ^M^ + ^ + ^n), 



rj=+l,-l 



(6) 



(7<X 



E jrJ dc i G> al { V + ^E + ^r)Ml{N n + l -- l -r ] ). 

T?=+l,-l ^ 



(7) 



The free generalized Kadanoff-Baym (FGKB) yj], [lOj ansatz for a two spin band system 
can be stated as, 



E a (p) + E al (p) 



(8) 



iG> a , (p, E, r, T) = -h(5^ - U (P, r, (E 



E ff (p) + ^(p). 



(9) 



The spin subband Wigner function is found by taking the energy integral of the G < , : 



/ <w /(p,R,T) = | dE(-i)G< CT ,(p,E;R,T). 



(10) 



Substituting the equations (@), © to equation (J2J together with the FGKB ansatz ©, 
© and using the expression (|Tn|) gives the desired set of equations. As a result, the spin up 



Wigner function for a Rashba effect RTD in weakly contact with a phonon heat bath can 
be written as 



~l Jteid P «co8[\z 1 (p,-p«)]( mmiz \^ ) + m . {z \ ^ )P^M\^T) 
-i- J d Zl d Pz2 cos[\ Zl { Pz -P^ ^ ¥) + mm{z \ ¥) )pJn(P\UP^z,T) 

+^ Jdz^sin^M-p ^^^ - m , ( , 1 +f) M/n(P|h^^,T) 
1 /" 1 

+^ / dz l dp z2 sin(-z l (p z - p z2 ))[E c (z - y) - £ c (z + -^)]/ n (p||,p z2 , z,T) 

+^Px(fn + /it) + ^Ps/(/n - /it) 

+ E E / - /t*(p + r, T))/ CTT (p, r, T) 

x.( ^ (P + q) ^- (P + q) - ^^(p) + tf^W, + I - i,) 
E E ^/rfq/T CT (P + q,r,T)(^ T -/ CTT (p,r,T)) 



<r=T,i »7=+i -i 



x ^ t( P ± q) + K(P + q) _ E AP ) + g ,(p) + + 1 + (n) 



The first five terms at right hand side of the above equation correspond to the drift terms. 
These terms are quite complicated comparing to the drift term in the usual single band 
Wigner function equation since the effective mass is taken to be position dependent. As a 
result, the effective mass becomes nonlocal. The sixth term gives the potential term and 
it is calculated self-consistently by coupling to Poisson equation. The seventh term is the 
subband mixing term. The last two terms correspond to the electron-phonon scattering. 
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The relaxation time approximation can be made for these. Similarly, the inter spin band 
polarization of the spin coherence becomes 



f d Zl d Pz2 cos[\ Zl { Pz -P«)\^ z { m . {z l _ f) + m * {z \ f) )pJn(P\UP^z,T) 

+^ Jdz^sin^M-p ^^^ - m . {z \ v )pifn(PhP«,z,T) 
If 1 

+ k / dZldPz2Sin (ft Zl ( Pz -Pz2))[E c {z - y) - £ c (^ + y)]/u(P||,P«2,^,r) (12) 

+|px(/u - /n) + «|pv(/tt - hi) 

+ E E Jfi / rfc l(^ - /t*(p + q, r, T'))/ (7 |(p, r, T) 
x , ( ^(P + 'q)+^(p + q) _ K(P)+^(P) + + 1 _ 1^ 

E E ^/rfq/Ta(P + q,r,T)(^-/ CTi (p,r,T)) 



<r=T,i »7=+l -1 



x ^ t(P + q) + K(P ± g) _ g^g) t g ,(p) + + 1 + ^ (13) 



Note that = ft,. Finally the spin down component of Wigner function is given as, 



k df n (p h p z ,z,t) 1 f . 1 1 1 v 2, / ^ 

* = 217 - P*)]( m . {z _ f) ~ m * (z + f) )p,fu(p\\,P«, z, T) 

X f X X X Q 

--] d Zl d Pz2 cos[- hZl ( Pz - Pz2 )]( m ^_ f) + ^^y^^/ubih^,-^) 

X /* X (9 X X 

-- J dz^cosi-zM-p^}- ^ _ + ^ g+ Wu(p|,,p^^T) 
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If 1 Z\ Z\ 

~Px(fn + /it) - i^Pyifn ~ /it) 

+ E E y rfc ^ - /^(p + ^ r ' T ))A(P' r > T ) 



xf( gl (p + q)+g,(p + q) _ ftW+MO + , nq)M , (JVq + 1 _ I,) 



- E E ^/ rf q/j<.(p + <J. r .r)(^] -/.j(p.r,r)) 

x ; ( ^pVq)^(p ± q) _ W*,(P) + „ q)Mq2(jVq + 1 + (14) 

The equations above can be further simplified by considering a constant effective mass 
and ignoring the scattering 



t.9f^(p\\,p z ,z,t) Pz d a a 

K df = -^^T + 2^Re[/ n ]-2-p 2/ Im[/ n ] 



+ h 

dRe[f n (p l] ,p z ,z,t)} p z d a 
h df = -m*Tz f " + h P ^- M 



If 1 Z\ Z\ 

- / dz 1 dp z2 sin{-z 1 {p z - p z2 ))[E c (z - y) - E c (z + -j-)]f n (p lb p z2 , z,T) 



If 1 Zi Z\ 

+ h / ^ 1 ^ 22 ^ n (^^ 1 fe~^ 2 ))[^ c ( 2; ~^")~ jBc ( 2; + y)]/Tl(^!l'^2,2;,T) 



dlm[f n (p h p z ,z,t)} p z d ,o> 

k df " ~^d~z f " + h Py{h ~ fll) 



If 1 Z\ Z\ 

+ h / ^ 1 ^ 22sin ^ Zl ^~ Pz2 ))[^ c ^~^")~ jEc ( 2; + y)]/iT(^ll'^2,^,T) 



If 1 Z\ Z\ 

+ h / ^ 1 ^ 22sin ^^ 1 ^~^ 2 ))[^ c ^~y)~ jEc ( 2; + y)]/^^!l' :P22 ' 2; ' T ) 
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where Re[/fjJ denotes the real part of fn and Im[/fjJ denotes the imaginary part of 

The particle density in each band is written in terms of the diagonal components of the 
Wigner function matrix as 

so that the spin up and spin down particle density are written as 



U| = j dp W dpz ^ p h Pz ' z ^ ( 16 ) 
H = ^3 J dp\\dp z f u (p\\,p z ,z). (17) 



The current density can be calculated using the following equation 



j M=E^/*^(p^). (18) 



Therefore the spin up and spin down current components can be written respectively as 



M z ) = v J dp \\ dp ^h(PhPz,z), (19) 

J i( z ) = J d P\\dPz^fu{p\\,Pz,z). (20) 

We derived the two-spin band Wigner function equations for a Rashba effect resonant 
tunneling diode. The effective mass became nonlocal. This provides the possibility of an 
accurate study of the relationship between the position dependent effective mass and the 
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spin splitting in the resonant tunneling structures since the interface effects on the effective 
mass are explicitly included in this model. Then we simplified the equations using a constant 
effective mass. This gives a set of equations that are simpler to solve numerically than the 
first set. The current and particle densities were derived. These equations are going to serve 
as the starting point for the Wigner function simulations of zero magnetic field resonant 
spin tunneling devices. 
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